Simplicity and scaling - size of a real polymer in three (or any) dimensions 
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We examine the scaling of the linear dimension of the system size of a real polymer solution 
at constant excess free energy and in two different spacial dimensionalities, d — do and d = di. 
Standard results for the functional form of the excess free energy lead to the conclusion that the 
scaling exponent i/(d) satisfies h'{do) ~ z^(cii) = dg^ — d'[^. Taking the critical dimensionality as a 
point of reference (i^(4) = 1/2) gives a scaling exponent — 1/4 + 1/d, in agreement with the 
accepted result for two- dimensions (i^(2) = 3/4) and the first term in the epsilon (d — 4) expansion. 
For the unsolved case of three dimensions it predicts i/(3) — 7/12. Several simplifying features of 
this result are pointed out. 



An amazing thing about polymers is that, although 
they are complex molecules, just a few parameters usu- 
ally describe their macroscopic properties [l[ . For exam- 
ple, the polymer size, Lp, is given by the scaling relation 
Lp ^ /q"^"*'''''- Here, the microscopic quantities are the 
number of monomers n, and a length ^o- The precise defi- 
nition of the latter we discuss later. The scaling exponent 
uld) generally depends on the spacial dimensionality d. 
If the monomers do not interact the chain is "ideal" and 
L ^ bn2 , independent of the dimensionality. Here b is the 
Kuhn length (the root mean square separation between 
adjacent monomers). More generally a polymer is "real" , 
meaning that the monomers do interact. Flory pointed 
out that this changes the scaling Using a very sim- 
ple model, he postulated that if the interaction between 
monomers is repulsive then, Lp ^ lQn'^Td{\ < d < 4) 
and Lp ^ Iqu^ {d > 4). So, for c? < 4 the real chain 
is expanded relative to the ideal chain. For d > 4 the 
chain scales in the same way as an ideal chain, so the 
interactions have no effect. In this respect Flory's model 
is correct. For two dimensions it predicts i^(2) — 3/4. 
Again, this is accepted as correct [3|. Nonetheless, it is 
regarded as inexact for the most important case - namely 
three dimensions. 

The reason for this is as follows. De Gennes famously 
showed the equivalence of the real polymer model to the 
n — limit of the n- vector model First, this shows 
that the real polymer problem has significance beyond 
the realms of polymer physics. A nice example is that 
modelling the universe as a real polymer resolves Olbers' 
paradox without assuming an expanding universe Q. 
Second, having identified this equivalence de Gennes 
could use Wilson's method [6] to show that near to the 
critical dimensionality {d = 4), Lp(e — > 0) ~ n'^ + TS, 
where e = 4 — d. In this limit (which we refer to as 
d — 3.99), the Flory result gives Lp{e 0) ^ lon'^ + 'h. 
Consequently, the prediction that for c? = 3, Lp ~ Ion's is 
approximate. Current numerical estimates are that ac- 



tually = 0.5877(±0.0007) '?]. So the exact value of 
this fundamental exponent, almost the polymer equiva- 
lent of 77, remains unknown. Here we use the functional 
form of the excess free energy of a polymer solution in 
different spacial dimensions to argue that in three dimen- 
sions 1^(3) = 7/12(= 0.5833). We also point out a number 
of simplifications regarding the physics of polymers and 
polymer solutions that would then follow. 

There are two microscopic lengths involved in the prob- 
lem, b and I, where I is the range of the monomer- 
monomer potential. For a given value of the ratio l/b, 
the size of the chain can be written in terms of either 
I or b. Here we take the former. Where proportionali- 
ties and similarities are used, from now on it is because 
we neglect constants dependent only on the dimension- 
ality. Consider the dimensionless excess free energy of 
a polymer solution, G*, where G* — G/kT with G the 
excess free energy, T the (constant) temperature, and k 
Boltzmann's constant. According to scaling theory Q 
this takes the form G* = Npfd{4>p), where fd is some 
dimensionally dependent function, Np is the number of 
polymers in the system and 4>p is the "polymer space 
fraction" 
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Here Lp is the linear dimension of the polymer in the 
dilute limit. In the scaling (n oo) limit this has the 
form 



Lp{d) - 



(2) 



Equally well, we can write the dimensionless excess free 
energy in terms of a dimensionally dependent function 
gd and the "polymer length fraction", as G* = Npgd{^p) 
where 



i/d-. 



(3) 
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Because this parameter involves lengths, which can be 
compared in different dimensionalities, it is more conve- 
nient for our purposes here. We now consider two sys- 
tems in two different dimensionalities do and di. The 
Kuhn length b in dimensionality do is proportional to the 
Kuhn length in dimensionality di, so there is no change 
in the relative size of the polymers in the two dimen- 
sionalities because of a varying ratio of Kuhn lengths. 
We now consider how the linear dimension of the sys- 
tem changes as we increase the number of the polymers 
and/or monomers, such that the length fraction remains 
unchanged. That is, such that G*{do)/G{di) depends on 
the dimensionality but nothing else. To do this requires 
that 



so the ratio of the system sizes scales as 



(4) 
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So the ratio of the system sizes in the two dimensionalities 
scales as 
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L{di) 
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The above form of the excess free energy also holds 
above the "overlap concentration" , $p > > 1 , where the 
polymers strongly overlap and the solution looks like a 
monomer soup. In this limit one expects that n is irrele- 
vant and only the monomer length fraction, $„, 



(7) 



is relevant. Here N,n{— nNp) is the total number of 
monomers in the system. In this limit the dimensionless 
excess free energy, in terms of a dimensionally depen- 
dent function hd, then takes the form G* — Nmhd{^m)- 
Far above the overlap concentration, but still in the 
semi-dilute regime ($p >> 1, $„i << 1) both func- 
tional forms of the excess free energy are valid, that is 
9d{'^p) ~ nhdi'^m)- This leads to generally accepted 
scaling of the free energy in this limit T] . Specifically, in 
terms of G* = A^p5d(*™n''^'*^"^/'*^ ), so for the above 
to be true requires G* jN^, ~ We should 

note, however, that for one dimension, where v(X) = 1, 
the dimensionless free energy is a function of only so 
both functional forms for the excess free energy cannot 
be satisfied simultaneously. This is reflected in the non- 
sensical prediction that in one dimension the free energy 
is infinite. 

Now we can again consider the transformation in di- 
mensionalities do and d\ (neither of which is now unity) 
described above. In terms of the monomer length frac- 
tion we have 
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In the semi-dilute limit, both [51 and 1101 are true, implying 
that 



v{do) - iy{di) ^ l/do - 1/di 



(11) 



If we take four dimensions to define do (that is do = 
4, 1^(4) = 1/2), then, because di can be any other dimen- 
sionality d, we have 



(12) 



For two dimensions this yields iy{2) = 3/4, the accepted 
result. For three dimensions it yields i/(3) — 7/12. Note 
that we could equally well take the two dimensional re- 
sult as a point of reference. In this case we would recover 
the four dimensional result. According to equation [T^] 
the excess free energy takes a particularly simple form 
above the overlap concentration, G* ~ Np^p, indepen- 
dent of dimensionality. We now consider a number of 
other simplifying features that result from equation 1121 

The monomer overlap concentration. Returning 
to the more usual space fractions, the thermodynamic 
properties of a polymer solution above the overlap con- 
centration are determined by the monomer concentra- 
tion, c ~ Nm/L'^. In terms of this variable the gen- 
eral form of the dimensionless excess chemical potential, 
^* = G* /Np, valid for all concentrations, is 



A* = /(c/ct) 
where c* is the monomer overlap concentration 
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According to equation [12] the monomer overlap concen- 
tration is c*~^ ^ (n^^'^b)'^. That is, the inverse monomer 
overlap concentration is proportional to a length (n^/'*6), 
independent of dimensionality, raised to the power di- 
mensionality. So, for polymer solutions the excess chem- 
ical potential, in terms of monomer concentration, takes 
the same form as that of a simple fiuid. For the latter, 
the length is the range of the potential (for example, the 
radius of the hard spheres in a hard sphere fluid). 

The dimensional dependence of the second 
virial coefficient. According to scaling theory the sec- 
ond virial coefficient of a polymer solution, B2, is pro- 
portion to the space occupied by the polymer. That is 
B2 cx Lp. From equation[T2l we then have B2/n (x b'^n'^^^ 
so the "second virial length" (^2/?^)^^'^ should be inde- 
pendent of dimensionality. In figure [1] we have plotted 
this quantity for the self avoiding random walk in two 
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case, we have 



FIG. 1: Second virial length (defined in the text), in two and 
three dimensions, as function of the number of steps of the 
seh avoiding walk. The dashed line is a fit to the asymgtotic 
two dimensional result. Data are taken from reference 



and three dimensions [7,]. The data for the two dimen- 
sionalities are in surprising good quantitative agreement. 
There is a small but statistically significant difference for 
smaller n, but for the largest value (n — 80000) the differ- 
ence is not significant. The reason that this is surprising 
is because our argument only predicts that the two func- 
tions are proportional. It appears from the simulation 
data, at least up to n = 80000, that the constant of pro- 
portionality is also independent of the dimensionality. If 
the result for the exponent given here is approximate the 
two curves shown in the figure must diverge for higher n. 
This could be tested numerically. 

Hyperbranched polymers. If an ideal chain with n 
monomers has m branch points, for large n and m its size 
is reduced by a factor oi m^^^ [ol], relative to the linear 
case, so 



r,l/4' 



(15) 



For a hyperbranched chain (m oc n) we therefore have 
Lp ^ n}/%. For all dimensionalities between two and 
four the hyperbranched ideal chain is an increasingly 
dense object. Its size increases more slowly with in- 
creasing mass than a space filling Euclidean object. 
This contrasts with the linear chain, which is rarefied 
in that its size increases more rapidly with increasing 
mass. Consequently, it is generally accepted that the 
effect of excluded volume on a hyperbranched chain is 
to change the sub-Euclidean scaling to Euclidean scal- 
ing, that is Lp [13] ■ On long length scales 
one expects that correlations can be neglected and that 
the real hyperbranched chain behaves like the equiva- 
lent hyperbranched ideal chain, but with an n dependent 
Kuhn length (the "uniform expansion" model). That is, 
Lp ~ ni/26„/mi/4, with 6„ = n''^<^^-'^/%. This being the 
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(16) 



For 



the 



hyperbranched chain we then have 
^Z"*, so substituting equation [T^ for the expo- 



nent gives the expected Euclidean scaling, Lp ~ n^^'^b. 

To examine the roles of the microscopic parameters, 
we now make use of the "two parameter model" JTi, . 
From dimensional analysis, the size of a single polymer 
takes the form 



Lrj 



(17) 



The unknown function f{z) is a function of a dimension- 
less (but dimcnsionally dependent) quantity z, 



(18) 



Now we have lifted the restriction that we are considering 
a system with a given value of the ratio l/b. 

The microscopic length. According to scaling the- 
ory Jj, for large z the size has a power- law dependence 
on n. That is, 



Lr 



^l-d7(d)^d7(d)^i+7(rf)(^) 
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with 7(d) = (2z/(d) - 1)/ (4 - d). So, the dimen- 
sional dependence of the scaling exponent also deter- 
mines the microscopic length determining the poly- 
mer size. For example, the Flory result gives = 
ld/(d+2)fj2/(d+2)^ Substituting equation [12] for the expo- 
nent yields 7((i) = l/(2d) meaning that 



Lp - l^b^n"'-'^'^ 



(20) 



There is a pleasing symmetry in the fact that the two 
intrinsic lengths in the problem, I and 6, enter with the 
same power. They play an equivalent role in determining 
the size of the polymer. The fact that is independent 
of dimensionality also means that, as for the ideal chain, 
the ratio of polymer sizes in two different dimensionalities 
is only a function of n. 

The excess free energy of dilute polymer so- 
lutions. Turning to the thermodynamics of low density 
polymer solutions (c/c* << 1), in the limit z ~> (where 
the polymer is hardly expanded) the excess free energy 
G* is 



G* - Nrad'' 



(21) 



Thus, the system is thermodynamically equivalent to a 
simple fiuid of monomers with size I. On the other hand, 
in the scaling {z — > 00) limit the excess free energy takes 
the form 



G* 



(22) 



In terms of the monomer concentration, using the two 
parameter expression for Lp, yields 

G* - iV™cz('^^('^)-i)Z'* (23) 

Introducing the the size of the polymer relative to the 
ideal chain size in the dilute limit, ao ~ z^'^'''\ equation!^ 
becomes 

i d-,{d)--i \ 

G* ^N„,ca)r^>f (24) 

Thus, in the scaling limit the excess free energy of a poly- 
mer solution, in terms of the monomer concentration, 
takes the same form as that for a simple fluid except that 
the interaction length is scaled by a factor related to the 
degree of expansion on the polymer. The dependence 
of this scaling on the dimensionality in turn depends on 
the dimensional dependence of the scaling exponent. For 
example, the Flory result gives G* ^ Nmcl'^/al, inde- 
pendent of dimensionality. On the other hand, using the 
expression for v{d) given by equation 1121 gives 

G* ^ N^c (^^^ (25) 

That is, the system is thermodynamically equivalent 
to a simple fluid, except that the effective size of the 
monomers is reduced proportionately to the degree of 
expansion of the chain, independent of dimensionality. 

As noted above, we cannot take one dimension as a 
point of reference because in this case the excess free en- 
ergy, above the overlap concentration, cannot be written 
in terms of both the polymer and monomer concentra- 
tion. It is unsurprising, therefore, that equation [T^] is 
incorrect for d = 1. It predicts iy{l) = 5/4. From the con- 
dition of fixed contour length, v^d) cannot exceed unity. 
One interpretation of this is that the expression predicts 
a lower critical dimensionality at d = 4/3 (the dimen- 
sionality for which the exponent is unity). We believe 
that there is tacit evidence for this from renormalization 
group theory (RG) calculations and simulations in 
non-integer dimensions [13] (see figure [2]). We use the 
word "tacit" here because ref. is, to our knowledge, 
only published as a pre-print and the author of ref. [ij] 
interpreted the low dimensionality behaviour as probably 
an error in the model because it disagreed with the Flory 
result (without explaining why the same model appar- 
ently worked for higher dimensionalities). 

Finally we should point out the reasons for believing 
that the result we obtain here is not exact. First, it is 
outside the range of current numerical and RG estimates 
for the exponent in three dimensions. However, the for- 
mer require extrapolation to the scaling limit and the 
latter the re-summation of divergent series. This makes 
reliably estimating the error difficult (when Nienhuis ar- 
gued that i'{2) — 3/4, this was outside the range of con- 
temporaneous numerical and RG values). Second, equa- 
tion [T2] agrees with the epsilon expansion to first order 




FIG. 2: Scaling exponent as a function of dimensionality 



but not to second order. Specifically, the epsilon ex- 
pansion gives z^(e) = 1/2[1 -I- e/8 -I- 15/256e^..], whereas 
equation [H gives i/(e) = 1/2[1 + e/8 + 1/326^..]. If the 
epsilon expansion is exact and unique to second order 
then, while equation [T^] is correct for d — 3.99 and d — 2, 
it is only a very good approximation for d = 3. Could 
the epsilon expansion to second order be inexact or non- 
unique? It seems heresy to even suggest this, but we 
note that there is an added degree of subtlety in calcu- 
lating the second order term in the epsilon expansion as 
compared to the first order term. Namely, a "magic in- 
teraction strength" is required for which corrections to 
scaling disappear [l5| . Further, the epsilon expansion by 
its nature requires the concept of non-integer spacial di- 
mensions. Are non-integer dimensions unique or just an 
interpolation between integer dimensions? In the analy- 
sis above we have treated d as a continuous variable, but 
this is not actually necessary. We could equally well have 
restricted ourselves to integer dimensionalities. The con- 
cept of non-integer dimensionality is not required here. 
To summarize, there are two possibilities. The epsilon 
expansion is exact and unique. All the above simplifica- 
tions apply in 4, 3.99 and 2 dimensions, but not quite in 
three. Alternatively, nature is kind and there is nothing 
special about three dimensions. 
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